Eccentric binary neutron star mergers 
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Neutron star binaries offer a rich phenomenology in terms of gravitational waves and merger rem- 
nants. However, most general relativistic studies have been performed for nearly circular binaries, 
with the exception of head-on collisions. We present the first numerical relativity investigation of 
mergers of eccentric equal-mass neutron-star binaries that probes the regime between head-on and 
circular. In addition to gravitational waves generated by the orbital motion, we find that the sig- 
nal also contains a strong component due to stellar oscillations (/-modes) induced by tidal forces, 
extending a classical result for Newtonian binaries. The merger can lead to rather massive disks on 
the order of 10% of the total initial mass. 



Introduction. Binary neutron star (NSNS) mergers are 
among the most promising sources of gravitational waves 
(GWs) for ground-based interferometers, as well as plau- 
sible candidates for the central enginge of short-gamma- 
ray-burst (sGRB). One of the main tools to study NSNS 
mergers is numerical relativity, which experienced a dra- 
matic development in the last 10 years achieving impor- 
tant results in NSNS [IH3] and recently also for black 
hole-neutron star (BHNS) simulations [H [5J 

In numerical relativity NSNS mergers have been al- 
most exclusively studied for circularized initial data, with 
the exception of head-on collisions [5j [7J (and see [5] for a 
nearly head-on Newtonian simulation with radiation re- 
action). Somewhat surprisingly, there are no numerical 
studies of orbits with eccentricity (discounting the resid- 
ual eccentricity [9] of quasi-circular data). One reason 
for this is an expectation that the dominant population 
of NSNS mergers in the universe results from evolution 
of primordial stellar binaries. For such systems, resid- 
ual eccentricity at the time of merger is expected to be 
low. For example, a recent investigation suggests that 
only between 0.2% and 2% of this class of mergers de- 
tectable by Advance LIGO /VIRGO will have eccentricity 
e > 0.01 (and e < 0.05) [TO]. With such small residual 
eccentricity quasi-circular templates will be able to de- 
tect the majority of these events [11] . and moreover tem- 
plates have been proposed for weakly eccentric binaries 
(e < 0.10) [T2] , However, recent studies have suggested 
there may be a significant population of compact object 
binaries formed via dynamical capture in dense stellar 
environments [TBI [T4"] . In [T2] , it is estimated that the 
rate of dynamical capture NSNS mergers may be large 
enough to account for a significant fraction of sGRBs, 
assuming NSNS mergers are their progenitors. Observa- 
tions of sGRBs have also suggested there are different 
progenitors for sGRBs that exhibit so called extended 
emission and those that do not |15) : an obvious specula- 
tion of the cause of the bimodality is primordial vs. dy- 
namical capture NSNS mergers. A large fraction of this 
latter class of mergers will occur with high eccentricity. 

Eccentric mergers can lead to a rich phenomenology in 
the GW and the merger remnant. The properties of the 
accretion disk formed as a product of the merger are so 



far poorly understood, even though there is the possibil- 
ity of massive disk production, M<j ~ 10% of the total 
initial binary mass. Recent results [16] for BHNS systems 
indicate a strong variability in properties of the merger 
remnant as a function of the eccentricity. The gravita- 
tional waveforms significantly differ from the chirp sig- 
nal of quasi-circular inspirals that feature slowly increas- 
ing amplitude. During eccentric orbits, each periastron 
passage leads to a burst of radiation, first studied us- 
ing Newtonian orbits together with leading order rela- 
tivistic expressions for radiation and evolution of orbital 
parameters [171 118) . and more recently with numerical 
simulations of BHBH and BHNS systems in full general 
relativity [16, 19-21]. An important Newtonian result 
concerning NSNS (or BHNS) is that eccentricity leads to 
tidal interactions that can excite oscillations of the stars, 
which in turn generate their own characteristic GW sig- 
nal [22] (see also [23, 24 ). In some cases the GWs can be 
dominated by these non-orbital contributions [22]. Nei- 
ther the orbital nor stellar GWs have been studied so far 
for eccentric NSNS orbits in general relativity. 

In this letter we present the first numerical relativity 
investigation of (highly) eccentric NSNS mergers. We 
consider an equal-mass binary at fixed initial separation 
and vary the initial eccentricity. Three models are dis- 
cussed that are representative of the different orbital phe- 
nomenology observed: direct plunge, merger after a close 
encounter, multiple encounters. We analyze the emitted 
GW, the dynamics in the matter and characterize the 
basic properties of the merger remnant. Dimensionless 
units c = G = Mq = 1 are used if not stated otherwise. 

Numerical Method. We performed numerical simula- 
tions in 3+1 numerical relativity solving the Einstein 
equations coupled to a perfect fluid matter model (no 
magnetic fields). We employed the BAM code [251125] , 
which we recently extended to general relativistic hydro- 
dynamics [27]. Referring to 26, 27J for details and further 
references, we solve the BSSN formulation in the moving 
puncture gauge coupled to matter in flux-conservative 
form. Metric and matter fields are discretizcd in space 
on 3D Cartesian meshes refined with the technique of 
moving boxes. Time integration is performed with the 
method-of-lines using a 3rd order Runge-Kutta scheme. 
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Derivatives of metric fields are approximated by fourth- 
order finite differences, while a high-resolution-shock- 
capturing scheme based on the local-Lax-Friedrich cen- 
tral scheme and the convex-essentially-non-oscillatory re- 
construction is adopted for the matter. Neutron star 
matter is modeled with a polytropic equation of state 
with adiabatic index T = 2. Gravitational waveforms 
are extracted at finite coordinate radius, r ~ 90M, by 
using the Newman-Penrose curvature scalar, \I/ 4 , where 
M = 2.8M is the gravitational mass of the system. 

Although the problem of posing initial data to the Ein- 
stein equations by solving the constraints is in principle 
well understood, there is no ready-made prescription for 
eccentric NSNS that are of interest here. We postpone 
the solution of the constraint equations and instead work 
with the following approximation. Initial data are set 
by superposing two boosted non-rotating stars with the 
same gravitational mass I.4M Q at apoastron. The two 
stars are located at coordinates (x±,y, z) = (±25M, 0, 0) 
and boosted with parameters £± = (0,±£ a ,0), respec- 
tively. The initial data are asymptotically flat but not 
conformally flat, and approximately irrotational. The 
maximum constraint violation (at the center of the star) 
decays as 1/d as expected for the relatively large initial 
proper separation d ~ 54. 7M. The constraint violation 
is at the level of the truncation error of the evolution 
scheme (see also the discussion in [16]). 

Orbital dynamics. We performed a series of simula- 
tions for different initial boosts (i.e. Newtonian apoastron 
velocities), £ a <E [0.01,0.05]. Newtonian circular orbits 
are obtained for £ y ~ 0.07. The grid configuration con- 
sisted of five refinement levels, maximum spatial resolu- 
tion h ~ 0.1 — 0.2M Q and a Courant-Friedrichs-Levy fac- 
tor of 0.25. For £ y < 0.020, within 1000M a merger and 
formation of a black hole occured. Selected models were 
continued for ~ 10ms after the merger in order to follow 
the accretion process onto the final black hole. During 
evolutions the largest violation of the rest-mass conser- 
vation occurs before collapse AMq/Mq ~ 0.01. Similarly 
the ADM mass was conserved up to 1%. The consistency 
of the results was assessed by convergence tests. 

We begin by reporting the orbital dynamics in our 
three models. See Fig. [T] for the star tracks as computed 
from the minimum of the lapse function |27j . 

In Model 1 (£ y = 0.010) the stars move almost directly 
towards each other, collide, then undergo prompt collapse 
to a black hole at t — r — 423M with dimensionless spin 
parameter obh ~ 0.58 and mass A/bh ~ 2.8Af Q . The 
resulting disk has negligible rest-mass, Md < 10 _7 Mo, 
which is at the level of the artificial atmosphere used for 
the numerical treatment of vacuum for the matter. 

In Model 2 (£ y = 0.020) the stars come in contact but 
survive the first encounter. During contact matter is ex- 
changed between the outer layers of the stars, see Fig. [2] 
which in turn gain rotation; the density weighted New- 
tonian vorticity is 10 times larger than an irrotational 
NSNS binary. During a second encounter they merge 
forming a black hole at t = 975M with initial obh ~ 0.79 
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FIG. 1: Star tracks for Model 1, 2, and 3. The tracks of the 
second star can be inferred from symmetry. 
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FIG. 2: Model 2 at t = 523M, shortly after the stars have 
touched and separated again. The rest-mass density (log- 
scale) and the three-velocity in the orbital plane are shown. 



and mass A/bh ~ 2.2M . In contrast to Model 1, there 
is a massive disk with Md ~ O.O8M0 ~ 0.15M©, mea- 
sured 100M after the formation of an apparent horizon; 
meanwhile accretion has increased the mass and spin to 
Mbh ~ 2.64M© and g,bh ~ 0.81. At the onset of merger, 
the separation remains nearly constant ~ 5M. This sug- 
gests a transition through a whirl regime as was observed 
in BBH systems P I UMH \ 125] . We stress, however, that 
the binary is in contact at this point and that the mat- 
ter interaction may lead to additional effects. After the 
merger we observe quasi-periodic oscillations in the accre- 
tion flow at a frequency isqpo ~ 0.5kHz in the vicinity 
of the black hole. 

In Model 3 (£ w = 0.022) multiple close passages with 
no matter exchange are observed. The two bodies do 
not merge during a simulation time of t = 5400 A/ . The 
orbits show a significant periastron precession, leading to 
an overall rotation of the quasi-elliptical orbit of almost 
90 degrees per encounter. However, zoom-whirl orbits 
are unstable, which suggests that by t = 5400M Model 3 
has not yet quite evolved to this regime. 

Waveforms. The orbital motion of each model results 
in a characteristic gravitational wave signal. Figs. [3] and 
[4] display the £ = m = 2 mode, ^22) as weu as the cor- 
responding instantaneous gravitational wave frequency, 
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FIG. 3: Real part of the r waveforms (central panel) and 
instantaneous gravitational wave frequency (bottom panel), 
M0J22, as a function of retarded time, (t — r)/M, for Models 1 
and 2. The top panel shows a zoom in on the shaded region 
in the central panel for Model 2. The horizontal line in the 
bottom panel marks the perturbative /-mode frequency. 
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FIG. 4: Same as the central and bottom panels of Fig. [3] but 
for Model 3. 



MW22, of the curvature scalar, which represents the main 
emission channel for the binary signal. 

Waveforms from Model 1 are characterized by a peak 
at the merger and the subsequent quasi-normal-mode 
(QNM) ringing of the final black hole. The GW fre- 
quency of the QNM is Mu>22 ~ 0.48, compatible with the 
fundamental QNM frequency of the Kerr black hole com- 
puted from the apparent horizon parameters, e.g. |29j . 
In fact, we find the value z^qnm ~ 5.9kHz, which differs 
by 6% from the value obtained from the GW. Waveforms 
from Model 2 show a burst at retarded time t — r^ 500M 
related to the orbital dynamics [IS], followed by a fea- 
ture between t - r ~ [800M - 900M] akin to a tran- 
sition through a whirl phase and the final merger sig- 
nal around t - r ~ 1000M . In this case the QNM fre- 
quency is higher, Mw 2 2 ~ 0.632 (z^qnm ~ 7.6kHz, 2% 



discrepancy). The waveforms from Model 3 exhibit sev- 
eral bursts corresponding to the periastron passages. 

The key feature is that between the bursts, (see e.g. for 
Model 2 at t - r ~ [550 - 750 M]). One can observe high 
frequency signals in several GW-modes, which are ab- 
sent in the black hole case, but qualitatively similar to 
BHNS eccentric mergers [IB]. This signature is progres- 
sively suppressed at lower eccentricities (larger £„) and 
not observed for £ y = 0.05. Referring to Fig. [4) the (2,2) 
GW frequency is dominated by several plateaus compati- 
ble with the /-mode frequency of the non-rotating star in 
isolation [5Drlg2] . Muf ~ 0.137 (y f ~ 1.586kHz), inter- 
rupted briefly by the close encounter bursts (the interup- 
tions at t ~ 1000M, 2000M disappear for larger extrac- 
tion radii). In Model 2, the signal is mainly in the (2, 2) 
multipole; the (2, 0) multipole (not shown) also contains 
the signature (weaker in amplitude by a factor five). In 
contrast, in Model 3 the amplitude of the (2, 2) multi- 
pole is smaller in amplitude by a factor of three than 
the (2,0) multipole, which is the main emission channel 
between the bursts (encounters). 

The natural interpretation is that these GWs are pro- 
duced by oscillations of the individual stars, which are 
tidally induced by the companion. Physically, the expec- 
tation is that the approximately head-on, zoom-part of 
the orbit is responsible for exciting axisymmetric m = 
modes, while in particular m = 2 modes should arise from 
the approximately circular, partial whirl near apoastron. 
The relative amplitude of the modes indicates the effec- 
tiveness of the different phases of the binary interaction 
to excite certain modes. Indeed the Newtonian analysis 
of Turner [32] indicates that a sufficiently close periastron 
passage (i.e. e < 1 for fixed apoastron) exerts a pulse- 
like tidal perturbation which excites the axisymmetric 
/-modes of the star, leading to a GW waveform domi- 
nated by the star oscillations rather then by the emission 
due to the orbital motion. Note that the phenomenon is 
different from a resonant tidal excitation, see e.g. [2"31 133] . 
which instead refer to the circular motion case. 

If these waves indeed correspond to stellar /-modes, 
they should be detectable as oscillations of the stel- 
lar matter. A way to identify the modes is to con- 
sider the projections of the rest-mass density of one of 
the stars onto the spherical harmonics, e.g. [31], pi m = 
I Yl m pd?x, and perform an analysis of the spectrum. 
Results are reported in Fig. [5] for Model 2 and 3, where 
the p-iiit) projections and their power spectral density 
(PSD) are shown. In both models we can clearly identify 
frequencies compatible with the linear /-mode (broad- 
ened by the signal length in Model 2), together with 
secondary peaks eventually due to nonlinear couplings. 
Similar results are obtained for the axisymmetric (m = 0) 
quadrupole mode, while in this case the signal is strongly 
modulated by the orbital motion. Overall the analysis in- 
dicates that both non-axisymmetric and axisymmetric /- 
modes are excited during the encounter. Note that, given 
an uncertainty on the frequencies, we can estimate an up- 
per limit on the rotation of the star. In model 3 there are 
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FIG. 5: Mode analysis. The panels display the I = m — 2 
spherical-harmonic projection of the rest-mass density of one 
of the stars (Model 2 top, Model 3 center), and its power 
spectral density (PSD, bottom), where Uf is marked. The 
width of the green shaded area, 73Hz, spans the perturbative 
values of the mode frequencies of models A0-A5 (right-left). 

more cycles leading to the most restrictive estimate for 
the maximum rotation rate (rotational energy/potential 
energy) of ~ 0.05, i.e. at most mildly rotating [30]. r- 



modes are thus not expected, but they may be present 
in the outer layers of Model 2; <?-modes are absent here 
by construction. 

Conclusion. Mergers of eccentric NSNS could be very 
interesting sources for third generation GW detectors, in 
particular given the (in some respect) surprisingly strong 
and clear signal from orbit-induced stellar oscillations 
(OISOs). If detected, this would pose constraints on the 
equation of state. The present results are likely to be 
highly dependent on the compactness of the stars and on 
the mass ratio, thus these cases should be studied as well. 
Furthermore, even larger disk masses are expected for the 
unequal mass case which could be interesting for sGRB 
astrophysics. More work is required in these directions 
together with a detailed understanding of the eccentric 
NSNS distribution in the Universe. 
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